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We consider a wave-function approach to the false vacuum decay with gravity and present a new 
method to calculate the tunneling amplitude under the WKB approximation. The result agrees with 
the one obtained by the Euclidean path-integral method, but gives a much clearer interpretation 
of an instanton (Euclidean bounce solution) that dominates the path integral. In particular, our 
method is fully capable of dealing with the case of a thick wall with the radius of the bubble 
comparable to the radius of the instanton, thus surpassing the path-integral method whose use 
can be justified only in the thin-wall and small bubble radius limit. The calculation is done by 
matching two WKB wave functions, one with the final state and another with the initial state, with 
the wave function in the region where the scale factor of the metric is sufficiently small compared 
with the inverse of the typical energy scale of the field potential at the tunneling. The relation 
of the boundary condition on our wave function for the false vacuum decay with Hartle-Hawking's 
no-boundary boundary condition and Vilenkin's tunneling boundary condition on the wave function 
of the universe is also discussed. 



In the previous literature, the false vacuum decay with gravity was investigated, in essence, by a naive extrapolation 
of well-established methods used for the Minkowski background jD . There are two methods to calculate the decay rate 
of false vacuum, one using the Euclidean path integral and the other constructing a WKB wave function. However, 
in both approaches such extrapolation is valid only for a tunneling potential V((j)) satisfying the thin- wall condition 
(we shall come back to this point below). 

The false vacuum decay in quantum field theory on the Minkowski background is analyzed in the semiclassical way 
(i.e., in the WKB approximation) as follows First, the Euclidean classical equation of motion for the field 4> (we 
consider a real scalar field throughout the paper) is solved by assuming the 0(4) symmetry on the field configuration 
with the boundary condition that the solution is regular at the origin and asymptotically approaches the false vacuum 
at infinite radius. There will be two such solutions; one trivially sitting at the false vacuum and the other that 
leaves the origin from a field value near the true vacuum and goes over the potential barrier. We call the former the 
false vacuum instanton and the latter the bounce instanton. We can take either the path-integral approach or the 
wave-function approach. 

In the path-integral approach, we sum up the contribution of instantons and fluctuations around them to obtain 
the tunneling amplitude or the decay rate of false vacuum. If we denote the Euclidean action of the false vacuum 
instanton by Sf and that of the bounce instanton by Sb, the decay rate is given by T ~ e~( Sa ~ SF \ In passing, we 
note that in the Minkowski background it is customary to set V(</>_f) = so that Sf = 0. But one can take a different 
choice, V((f>F) ^ 0, in which case both Sf and Sb would diverge but the difference is well-defined and independent 
of the value of V(4>f)- 

In the wave-function approach, we construct a one-parameter family of 3-dimensional field configurations from the 
everywhere-false- vacuum configuration to the configuration with the true vacuum bubble about to expand in the false 
vacuum sea (the critical bubble configuration). Then, exponentiating the action of these configurations we obtain the 
WKB wave function(al) describing the tunneling of the field from the false vacuum state to the state of the critical 
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bubble H . Then the tunneling amplitude is obtained by taking the square of the ratio of the final state wave function 
to the initial state wave function, which is also given by e~( SB ~ Sp \ 

Since these two different approaches give the same result, and since there is one-to-one correspondence between a 
classical solution and a WKB wave function, we may justify the interpretation of the bounce instanton as describing the 
tunneling process and its analytic continuation to the Lorentzian space-time as describing the classical configuration 
of the field after false vacuum decay. 

A trouble arises when one attempts to extrapolate the above procedure to the case of false vacuum decay with 
gravity. First, the scalar field equation is now coupled with the Friedmann equation that governs the scale factor of 
the 0(4) symmetric metric. Second, and most problematically the 4-dimensional Euclidean space becomes inevitably 
compact with topology of S 4 when gravity is taken into account (we assume V(</>f) > where <pF is the field value 
at the false vacuum and a potential of the shape depicted in Fig. Ill) , i.e., not H 4 as in the case of flat space-time. Of 
course a bounce solution still exists in the presence of gravity, first obtained by Coleman and DeLuccia [jjj. However, 
there will be no bounce solution that admits a 3-dimensional slice on which the field is asymptotically at false vacuum 
everywhere; an instanton leaving, say, the north pole of S 4 from the true vacuum side does not reach the false vacuum 
when arriving even at the south pole. Furthermore, the standard dilute gas approximation to obtain the decay rate in 
the path-integral approach, in which the 4- volume occupied by a vacuum bubble is assumed to be negligible compared 
with the whole Euclidean 4-volume, fails. 

Nevertheless, there is one special case in which the extrapolation can be almost justified. It is the thin-wall limit 
|^| . In this case, the field </), which is sitting on the true vacuum side of the barrier at the north pole, varies abruptly 
at the position of the bubble wall sharply located in the northern hemisphere and cf> — cf>F everywhere outside the 
bubble. Then, the solution admits a maximal 3-surface on which the extrinsic curvature of the 3-surface vanishes 
and <f> = <f>F everywhere. Through this surface, we can analytically continue to the Lorentzian solution that describes 
the false vacuum state. The critical bubble configuration that describes the moment of bubble nucleation is also a 
maximal 3-surface, and analytic continuation through the surface gives the space-time with an expanding bubble. 
Thus we obtain a one-parameter family of 3-dimensional configurations that interpolates between the false vacuum 
and critical bubble configurations. 

However, only a very restricted class of potentials admits bounce instantons in the thin-wall limit. Namely, the 
ratio of the mass scale m of the curvature at the top of the potential barrier, m 2 = —d 2 V/d(j) 2 ((f)top), to the typical 
mass scale M of the potential energy density must be very large (rigorously speaking, it must be infinitely large in the 
exact thin- wall limit). In other words, the barrier must be extremely sharply peaked. We note that this is somewhat 
different from the definition of the thin- wall limit originally discussed by Coleman and DeLuccia Q| . In our definition, 
since the Euclidean 4-volume is proportional to (GM 4 )~ 2 while the 4-volume of a bubble is 0(m~ 4 ), the thin-wall 
limit implies small bubble radius, provided M is smaller than the Planck scale, i.e., GM 2 < 1. Field potentials with 
such a feature are apparently not general. However, previous discussions on the false vacuum decay with gravity have 
been relied on the picture of a bounce instanton having the feature of the thin-wall limit. For example, it is widely 
believed that the formula T ~ e -( s B-s F ) j g g^jrj valid even when the bounce solution does not have the thin-wall 
feature at all. Furthermore, the interpretation of such a bounce as describing the tunneling process and describing 
the classical evolution from the critical bubble state has been adopted without serious considerations. 
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FIG. 1. A schematic picture of the potential V(4>) of a scalar field discussed in the text. 
In this paper, we propose a new method to calculate the wave function describing false vacuum decay with gravity, 
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which is not restricted to the thin- wall limit. Our method gives a picture of tunneling substantially different from 
that in the Minkowski background, but nevertheless our wave function interpolates between the false vacuum state, 
which is described by the trivial instanton sitting at false vacuum, and the state at the bubble nucleation, which is a 
maximal 3-slice containing the critical bubble of the bounce solution. 

The paper is organized as follows. In Section II, we first review the mini-superspace Wheeler-DeWitt equation for 
0(4) symmetric configurations and the WKB approximation. Then we describe our method. We construct a one- 
parameter family of spatial configurations that interpolates between the false vacuum state and the critical bubble 
state. As we are interested in the WKB approximation, we consider a classical path connecting these two states. As 
pointed out in the above, no 'single' classical solution admits such a path. However, a unique feature of false vacuum 
decay with gravity is that one can match the 'two' instantons (the false vacuum instanton and the bounce instanton) 
smoothly across the 'south pole' of each instanton. In section III, by solving the Wheeler-DeWitt equation we explicitly 
perform this matching and construct the wave function that contains both the false vacuum and critical bubble states, 
hence describes the false vacuum decay with gravity for a wide class of tunneling potentials. It. is then straightforward 
to calculate the tunneling amplitude. The result recovers the Minkowski result if we take the zero-gravity limit G — > 0, 
and also agrees with the decay rate obtained in the path-integral method with a naive extrapolation of the formula, 
r ~ & -{Sb-S f ) _ Q ur resu it supports the standard interpretation that the bounce solution, despite the fact it does 
not contain the false vacuum configuration at all, does describe the classical evolution after false vacuum decay by 
analytic continuation through the critical bubble configuration. Section IV is devoted to discussions. A covariant 
formulation of the WKB wave function for tunneling is recapitulated in Appendix A. 

In the rest of the paper, we keep h and G explicit to clarify the WKB order as well as the effect of gravity. 



II. FORMULATION 



In this section, we describe our method to construct a relevant classical path in the Euclidean regime that determines 
the WKB tunneling wave function. 



A. Instantons and WKB approximation 



We consider the action, 



5' 



d i x^/g 



n 



24ir 2 n 



(1) 



where k := 2G/3ir, and the potential V(<fi) as illustrated in Fig. [j]. We confine our attention to WKB wave functions 
that are described by classical solutions having 0(4) symmetry. Hence we consider the metric in the form, 



ds 2 



-N 2 (t)dt 2 +a 2 {t)dnf 3) 



(2) 



in which df2 2 3) is the metric on a unit 3-sphere, and (f> = 4>{t). In this case the action in the first-order form becomes 



S 



dtL ; L = P n a - 



~2a 



-P 2 



(2T7 2 n)a 2 



P 2 



Then the Wheeler-DeWitt equation is written down as 

d 2 l2 p 3 1 

h -— 



K 

2a 



da 2 



a da (2tt 2 k) a 



2 o- or 



* = 0, 



(3) 



(4) 



where 



Q(a,(j)) = 1- (47r 2 K)a 2 V(^), 



(5) 



and p is an arbitrary constant that determines the operator-ordering; p = 1 corresponds to making the ordering 
covariant with respect to the superspace metric (see Appendix A). 

It should be noted, however, that it is not always necessary to take the t =const. 3-geometry as the argument of 
the wave function. In fact, the critical bubble configuration that describes the moment of bubble nucleation does not 
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respect the 0(4) symmetry. Nevertheless, we shall see below that the tunneling amplitude can be calculated solely 
with the knowledge of the wave function for 0(4) symmetric 3-geometries. 

The construction of a WKB wave function describing multi-dimensional quantum tunneling was discussed much 
in the literature || and reformulated in a covariant manner in 0]. Given a one-parameter family of configurations 
that satisfies the Euclidean equations of motion, the WKB wave function along the family can be obtained with this 
method, as recapitulated in Appendix A. Following this method, we first express the wave function in the form, 

* = exp[-W (a, <t>)/n - Wiifl, <f>) +■■■}■ (6) 
Then at the lowest WKB order, we have 

which is, of course, free from the operator-ordering ambiguity. 
Introducing a parameter r such that 

Mrl = _kOWo and #(t) ^ 1 8 Wp 

dr ' a da dr 2-K 2 a 3 50 

we find 

^-W (a,4>)=--Q(a,<fi). (9) 

The Euclidean equations of motion for 0(4) symmetric configurations, (o(r), </>(t)), with the choice of the lapse N = 1 
are obtained from Eqs. (0), (|j) and (||) as 

a 2 - (2n 2 n)a 2 ^ 2 - Q(a, 0) = , 

+ 3-0= — , 

a dip 

ad + 2d 2 + 1 = 3Q(a,4>). (10) 

The first equation is the (Euclidean) Friedmann equation that corresponds to the energy integral of the remaining 
two equations. The regularity of the solution requires = at a = 0. 

Equation (|J) readily gives the relative amplitudes of the leading-order WKB wave function at two different config- 
urations at t = T\ and T2, 



*(7s)=exp[-Wb(7a;ri)]*(ri); W (t 2 ;t 1 ) = — / dr' a(r')Q (a(r'),0V)) ■ (H) 



T2 



As is well-known, Wo(t2\t\) is just the Euclidean action integral of the system from the configuration at t\ to that 
at r 2 . 

The solutions relevant to our discussion are the trivial solution sitting at the false vacuum and the bounce instanton. 
The trivial instanton is explicitly given by 

(a FV »,0 FV ») = (H^MHft),^) ; t&[0,tt/H f }, (12) 



where Hp — kV(4>f)- It represents a Euclidean 4-sphere of radius Hp 1 and we call it the false vacuum (FV) 
instanton hereafter. The bounce instanton is known as the Coleman-DeLuccia (CD) instanton We denote it 

by (a (t), c - d (t)) with r £ [ti,T2\. It also has the topology of 5* 4 but the field varies monotonically over the 
potential barrier as the scale factor a varies from zero to a maximum, then to zero. For later convenience, we set 
(a CD (Ti), c - d (ti)) = (0, 0i) where 0i is on the true vacuum side of the barrier, and {a c D fa) , (j) C D fa)) = (0,02) 
where 02 is on the false vacuum side of the barrier. See Fig. |l|. The FV and CD instantons embedded in H 5 are 
schematically shown in Fig. || with 2-dimensions suppressed. 
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FIG. 2. Schematic pictures of the false vacuum (FV) instanton and the Coleman-DeLuccia (CD) instanton with 2-dimensions 
suppressed. The arrow indicates the direction of the increase in the parameter £ defined in the text. 

B. One-parameter family of tunneling configurations 

As we noted in Introduction, the CD instanton itself does not admit a maximal 3-surface on which <p = far 
everywhere, unless the potential satisfies the thin-wall condition. Hence the CD instanton alone cannot describe the 
tunneling wave function. It is then natural to suppose that the FV instanton plays a role as well. In fact, there is 
already a hint in the formula for the decay rate, T ~ e -(SB-SF)/h^ v& \[^ m th e thin-wall limit, where the action of 
the FV instanton comes into play. 

We therefore consider a possibility to construct a one-parameter family of spatial configurations that connect 
the false vacuum state and the critical bubble state by matching these two instantons somehow. To look for this 
possibility, we have to keep in mind that the WKB approximation requires us that such a path should satisfy the 
classical equations of motion almost everywhere. But this is the crucial point; it is so not everywhere but almost 
everywhere. In any calculation of a wave function under the WKB approximation, there can be configurations of 
measure zero that violate the WKB condition, such as turning points of a classical solution, but their existence does 
not invalidate the calculation if one appropriately performs the matching, e.g., by the method of asymptotic matching. 
Then we realize that the two instantons may be actually matched through the point a = where the 3-geometry 
ceases to exist, hence can be regarded as a turning point of the classical solutions. 

With the above considerations in mind, we combine and re-organize the 0(4) symmetric families of configurations of 
the two instantons discussed in the previous subsection to a one-parameter family of configurations which adequately 
describes the false vacuum decay. We denote this one-parameter family by (/iy (a;*; £), ^(x 1 ; £)) where £ is a non- 
dimensional parameter. The parameter £ is assumed to run through a range [0,£/] with the initial (£ = 0) and final 
(£ = £/) configurations given by the the maximal 3-slice of the FV instanton and the critical bubble configuration of 
the CD instanton, respectively. 

The initial configuration, denoted by E$ in Fig. ||, is given by Hpr = n/2 of Eq. (|l2|). We thus set £ = HpT — n/2 
until the scale factor a vanishes, 

{hijix^&fa'iQ) = (i/ F 2 sin 2 (£ + |)^ 3)i F ) ; £e [0,tt/2). (13) 

We match a — of the above to the point (a CD , <fi CD ) = (0, fa) of the CD instanton, where fa is the field value 
at the 'south pole', see Figs, [j] and ||. Note that the values of <f> are different at these two points where the WKB 
breaks down. However, as we shall see in the next section, this will not be an obstacle since the wave function will be 
independent of <p in the vicinity of a = 0. 

For £ > 7r/2, we take the 0(4) symmetric configurations of the CD instanton up to the slice S m shown in Fig. 
where E m is the 0(4) symmetric slice at a = where the 3- volume is maximum, which we denote by (hij, <f>) — 
(H m 2 dfl 2 {3) ,cb m ). Thus 

{k j (x i ;0,^ i -,0) = (a CD (^ts),<P CD (0)-, £e(~£ m ], (14) 

where 

(a CD (7r/2), <t> CD {n/2)) = (0, fa), (a CD (U),<t> CD (U) - (H m \<t> m ). (15) 
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For £ > £ TO , we slice the CD instanton in such a way that all of the configurations contain the common 2-sphere 
the slice S m and the final slice £/ of the critical bubble configuration intersect. These slices do not respect the 0(4) 
symmetry but are analogous to the static slices of the de Sitter space of radius H~ x : 

dr 2 

As mentioned above, we label the final configuration by £ = £/. Thus 

(hij(£f), </>(£/)) = ( tne critical bubble 3-geometry). (17) 

As it is clear from the above construction, the lapse function N vanishes at £ m , and the geometry is degenerate 
there. However, this will not be a problem since we have started from the Wheeler-DeWitt equation (^), in which 
the dynamical variables are {hij{x l ) 1 4>{x 1 )) of a 3-geometry. Hence, although the lapse function plays a role when a 
Euclidean classical solution is considered, the vanishing of it is irrelevant to our discussion (in this connection, see |Q). 
Furthermore, as the relative magnitude of the WKB wave function at two different configuration is simply determined 
by the action integral between them, and it is independent of deformation of a path between them (as long as the 
path lies on the space of the classical solution that dominates the wave function) , the contribution from the upper- left 
quadrant of the CD instanton in Fig. || swept both by the configurations between £ 6 (tt/2, £ m ] and £ € [£ m , £/] cancel 
each other and only the action of the lower hemisphere determines the relative magnitude of the wave function at 
£ = 7r/2 and £ = £/. Thus, although the values of the wave function for configurations in the range £ € (£m,£/) are 
difficult to obtain, they are not needed in the calculation of the tunneling amplitude. 



III. TUNNELING WAVE FUNCTION 



Let us now construct the tunneling wave function. Since we are interested in the false vacuum decay, we want the 
wave function to describe the expanding universe after tunneling. Thus the appropriate boundary condition for the 
wave function is the tunneling boundary condition, which demands the wave function describing the classical universe 
after tunneling to have a positive eigenvalue for the momentum operator P a = —ifid/da 

In accordance with the construction of the classical path given in the previous section, we assume the universe after 
tunneling is described by analytic continuation of the CD instanton through the critical bubble configuration. We 
choose the moment of nucleation (£ = £/) as t = t = 0. Thus we set 



= ^exp 



t(0 



dt'L CD {t')+i- 



(18) 



where TV is a normalization constant and L CD (t) is the Lagrangian of the Lorentzian CD solution. The phase ir/4 is 
inserted for convenience, and we have neglected the prefactor arising from the next WKB order. Then, the standard 
WKB connection formula gives the under-barrier wave function for £ G (7r/2,£/] as 



dr'L c E D {r') 







+ l -N exp 





\ f dr'L% D (r') 

h Mi) 



(19) 



where L^ D (r) is the Euclidean Lagrangian of the CD solution. 

For £ e (tt/2, £ m ] where the 3-geometries are 0(4) symmetric, we have from Eq. (O 



(0 



dT'L% u {r')= Wb(£/;7r/2) - W (^tt/2) 



Sb 



Wo(£;7r/2), 



(20) 



where Sb is the total action of the CD bounce instanton and 



Wo(f;7r/2)=-- / dr'a CD {r l )Q(a CD {r'),cp CD {T 1 )) 

K Jt(tt/2) 

=■■ w Q CD (0, 



(21) 



G 



where a C£> (r(7r/2)) = 0. Hence 

*(tt/2 < £ < £ m ) = ,4 exp 
where the coefficients A and B are given by 



+ Bexp 
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Likewise, the WKB wave function for £ € (0, 7r/2) is expressed as 

I. 



*(0 < £ < tt/2) = Cexp 

where C and D are constants and 

W FV (0 :=Wb(7r/2;0 



I? exp 



(22) 



(23) 



(24) 



Tv/Hjr 



(f) 



(25) 



where a Fy (r) is given by Eq. (jl^) . 

We consider the matching of the above wave functions in a region a ~ 0. Noting that V(<p) > V(0p) in the vicinity 
of £ = 7r/2, it is convenient to introduce a re-scaled (non-dimensional) scale factor, 



a := e- 



where the parameter e is defined by 



e := 2tt* 



>n 6 V{(j> F ) 



(26) 



(27) 



and £ p i := ^/Jik — ^j2hG/3ir is the Planck length and m p i :— \fhj~K — ^/3nh/2G the Planck mass. We assume the 
potential energy scale is mu ch below the Plank scale, e <C 1. 

In the region a <C 1 / (2n t/kVJ^f)) , which corresponds to a <C 1, Q(a, 4>) is independent of </>; Q « 1, but the WKB 
approximation is valid as long asa>e. Hence, Wq D ((,) and W<f v {£,) in this region reduce to the same form, 



The prefactor due to the first-order WKB correction is calculated to be (see Appendix A) 

p+i 

exp(— W\) = a 2 x const.. 



(28) 



(29) 



Although we are not interested in this first order correction, we recover it in the region e < o < 1 to show explicitly 
that the matching we perform below is independent of the operator-ordering. Thus the WKB wave function in the 
region e <C a <C 1 becomes 



and 



> 7r/2) = Ad ~ exp 



< tt/2) = C~aT— exp 



a 2 (0 

2e 2 



a 2 (0 



2r 2 



D ~_2±l 

±>a 2 exp 



p+i 

-I- Z?a 2 exp 



2c 2 



2e 2 



(30) 



(31) 



To match these two wave functions, we directly solve the Wheeler-DeWitt equation in the region < a -C 1 and 
compare the coefficients in the region e <C a <C 1. See Fig. |^. In terms of a, the Wheeler-DeWitt equation (Q) reduces 
to 



7 



p d a 



7,2 " 



da 2 



* = o. 



(32) 



The general solution is 



h 2 a 2 



(33) 



where v = (1 — p)/4, and F are constants, and /^(z) and K v (z) are the modified Bessel functions. The asymptotic 
form of this wave function at a 3> e is given by 



*(M) = 



eisa 2 exp 



2e 2 



exp 



2e 2 



2e 2 

(l + 0(— )). 

ar 



(34) 
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FIG. 3. A picture of the superspace (a, </>) near the matching region. The thick lines are the one-parameter family of 
configurations (foij(£)> $(£)) m ^ ne range £ £ [0,7r/2) and £ € (7r/2,£/]. The arrow indicates the direction of increase in £ along 
the family. The WKB approximation is valid in the hatched region and the wave function ( |34| ) is valid in the grayed region. 
The overlapping region is where the matching is performed. 



Since Eqs. ( |30| ) and (|31j) have exactly the same a dependence, the comparison of them with ( p4| ) readily gives 

A = eE = C , B = e(EeS u+ ^ m + F) = D , (35) 



independently of the choice of the operator-ordering. 

Plugging this result into Eq. (p4|) and noting Eq. (p3|), we find the wave function for £ G (0,7r/2) as 



*(0 < £ < tt/2) = Ne SB / 2h exp 



+ -Ne~ SB/2h exp 



(36) 



Noting the fact Wq (0) = Sf/2, we connect this to the Lorentzian region £ < of the false vacuum to obtain 











N exp 




jexp 


~ 1 

ft Jo 



-N exp 



2h 



(Sb — Sf) 



exp 



i 

ft Jo 



dt'L FV (t') 



exp 



dt'L FV (t" 



i 

n 

exp 



i 

ft Jo 



dt'L FV {t') 
*(«) 



dt'L FV (t') 



(37) 



Since Sb > Sf, the first two terms in the curly brackets dominate. Thus, the tunneling amplitude T of the false 
vacuum decay is found to be 



l*(0)| 2 



exp 



1 



(S B - S F ) 



(38) 



This agrees with the result in the thin- wall limit. The overall behavior of the wave function is shown schematically in 
Fig.^J, in which the normalization constant N is chosen to be exp[— (Sb — Sf)/2K\ so that the amplitude of the wave 
function at the false vacuum equals unity. 
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FIG. 4. A schematic behavior of the tunneling wave function. The point £ = n/2 is where the scale factor a vanishes, and 
towards both right and left directions from there the 3-volume increases. The dotted curve indicates the effective potential 
barrier. 



IV. DISCUSSION AND CONCLUSION 

We have presented a new method to calculate the tunneling wave function that describes the false vacuum decay with 
gravity. In this method, the tunneling wave function is constructed by matching the false vacuum (FV) instanton and 
the Coleman-DeLuccia (CD) instanton through the point at which the scale factor vanishes. We found the resulting 
tunneling amplitude agrees with the naive extrapolation of the formula whose validity had been justified only in the 
thin-wall limit or in the flat background. Our result is a strong support for the standard interpretation that analytic 
continuation through the critical bubble configuration of the CD instanton describes the universe after false vacuum 
decay. 

In our picture, the false vacuum decay with gravity becomes more or less quantum cosmological. For example, as 
can be seen from Fig. |], interpreting the wave function as describing the actual tunneling path, the tunneling under 
consideration is quite similar to a case of standard quantum tunneling with a resonant state inside the barrier, except 
for the sign of the Euclidean actions Sf and Sb which are negative. The difference is that the resonant state here 
is not a state of something but the state of 'nothing'. This gives us a very interesting picture that the false vacuum 
decay is first proceeded by pumping up the amplitude in the state of nothing and subsequently a universe with an 
expanding bubble is created from nothing. 

A similar idea of joining the two instantons was previously proposed by Bousso and Chamblin in the context of 
path-integral approach |||. But the principles behind his approach and ours seem rather different. Our method is 
completely within the scope of the standard WKB approximation, while they had to introduce a Planck size wormholc 
and reverse the sign the 4-volume after traversing the wormhole in an ad hoc way. Nevertheless, it is of interest to 
see if these two methods correspond to each other at certain level of approximations. 

Recently, Rubakov and Sibiryakov proposed another method to deal with the false vacuum decay with gravity ||, 
by considering complex paths and adding a constraint to the action to realize the initial false vacuum state. It is, 
however, not clear how the complex path they considered joins the false vacuum and the critical bubble configuration, 
particularly because of the constraint that modifies the equations of motion. 

We have demanded the wave function to have only the component describing an expanding universe after false 
vacuum decay In this respect, our boundary condition is very similar to the tunneling boundary condition proposed 
by Vilenkin |q{ in the context of quantum cosmology. However, in Vilcnkin's picture, only the wave function for an 
expanding universe appears in the Lorentzian region of the superspace. This means there is a steady flux coming 
out from the Euclidean boundaries of the superspace. On the other hand, our wave function is outgoing only for the 
universe after false vacuum decay. Indeed, along the line <fi — 4>f in the superspace, where the universe is in the false 
vacuum, our wave function has both expanding and contracting components. See Fig. |^. The latter feature is similar 
to the one of the Hartle-Hawking wave function . These similarities of our wave function with both the Vilenkin 
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and Hartle-Hawking wave functions can be seen also in Fig. [|: If we consider only the right-hand side (£ > 7r/2) of 
the figure, it looks exactly like the cosmological wave function satisfying the Vilenkin boundary condition. In fact, the 
amplitude at a — (£ = 7r/2) is dominated by the component that decreases exponentially as a increases, as is the 
case of the Vilenkin wave function. On the other hand, if we consider only the left-hand side (£ < ir/2) of the figure, 
our wave function can be approximately regarded as a Hartle-Hawking wave function. The exponentially decreasing 
component becomes totally negligible when it appears in the Lorentzian region, hence it is an irrelevant component, 
while the exponentially increasing component dominates the wave function in the Lorentzian region, which consists 
of both outgoing and ingoing components with equal weight. This is the characteristic Hartle-Hawking feature. 




FIG. 5. Schematic picture of the causal structure of superspace and the one-parameter family of configurations used 
to construct the wave function. The Euclidean region is shaded. The lines indicated by FV and CD correspond to the 
configurations of the false vacuum and Coleman-DeLuccia solutions, respectively. The up or down arrows along the lines 
indicate the expanding or contracting components contained in the wave function. 



We have considered only the leading order behavior of the WKB wave function. We need to analyze the next order 
to determine the quantum state of fluctuations after tunneling. It seems, however, that there is no reason to expect 
the result should differ from the ones obtained previously under the assumption of the Euclidean vacuum associated 
with the CD instanton. This is because all the properties of the quantum fluctuations are determined by the Euclidean 
structure of a single classical solution, which is the CD instanton, as far as we are interested only in the state after 
false vacuum decay: The quantum fluctuations around the CD instanton is completely decoupled from those around 
the FV instanton. Then once we focus only on the CD instanton, our wave function satisfies the outgoing boundary 
condition, and it is known to give the Euclidean vacuum state @. This reassures the validity of previous calculations 
of quantum fluctuations in the one-bubble open inflation scenario Jl^ |. 
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APPENDIX A: COVARIANT FORMULATION OF MULTIDIMENSIONAL TUNNELING 



Here we recapitulate the covariant formulation of the WKB wave function for multidimensional quantum tunneling 
developed in [Q, slightly adapted to the Wheeler-DeWitt equation. 
We consider the Hamiltonian in the form, 



H = -h 2 g a ' 3 (q)-^-V a f(q)S7 p + U(q) 

fw 



(Al) 



where {q a } are the coordinates in the configuration space of the dynamical variables, or superspace, and g a p is the 
superspace metric. Expanding the wave function as 



= exp 



1 



(W + hWi 



(A2) 
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and inserting it into the Wheeler-DeWitt equation H^> = 0, we obtain the lowest-order WKB equation, 

- g aP V a W o V W Q + U(q) = 0, 

and the first-order WKB equation, 



2 ff /(?) 



Introducing a parameter t such that 



d -^ = g^w , 



we obtain from Eq. (A3), 



W = / dTU{q{r)), 



and from Eq. (A4), 



Wi = - In 



det 



g\f(q) 



const., 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



where {A Q } = {r, A"} and {A™} labels the different orbits of the congruence satisfying Eq. (A5) on the superspace. 
For our mini-superspace system (Q), the choice {q a } = {a, <p} leads to 



.'/.. < I n oJ„3 ) . f(q) = a p \ 



The congruence containing both FV and CD instantons is expressed in the vicinity of a = as 

a(r,0 o ) =r + 0(r 3 ), 0(r, O ) = 0o + 0(r 2 ). (A9) 



Hence a conv enient choice of the label A of the orbits is the value of <f> at r = 0, <f>o. Substituting Eqs. (A8) and ( |A9| ) 
into Eq. JA7|), we find 



W\ = In a 2 + const. 



Equations ( |A6| ) and (A7) give the WKB wave function to the first order, 

* = Na- 2 ^ exp[- y dr [/(a(r), 0(r))], 

with J7 = —(a/ ' K)Q(a, </>). 



(A10) 



(All) 
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